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Abstract 

Motivated by recent experimental progress to read out quantum bits implemented in supercon- 
ducting circuits via the phenomenon of dynamical bifurcation, transitions between steady orbits 
in a driven anharmonic oscillator, the Duffing oscillator, are analyzed. In the regime of weak dis- 
sipation a consistent master equation in the semiclassical limit is derived to capture the intimate 
relation between finite tunneling and reflection and bath induced quantum fluctuations. From the 
corresponding steady state distributions analytical expressions for the switching probabilities are 
obtained. It is shown that a reduction of the transition rate due to finite reflection at the phase- 
space barrier is overcompensated by an increase due to environmental quantum fluctuations that 
are specific for diffusion processes over dynamical barriers. Moreover, it is revealed that close to the 
bifurcation threshold the escape dynamics enters an overdamped domain such that the quantum 
mechanical energy scale associated with friction even exceeds the thermal energy scale. 
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I. INTRODUCTION 



The prospect to tailor devices for quantum information processing has stimulated major 
experimental research in the last years. Several different technologies have been explored 
to assess the possibility to realize quantum bits, like ion traps liquid state magnetic 
resonance, linear optics, electrons in liquid helium and superconducting Josephson Junction 
(J J) devices (2-4]. Particularly for the latter ones the insulation of the structure in which the 
quantum bit is implemented, e.g. a Cooper pair box, from its surrounding is of substantial 
relevance. This issue also includes the readout device of the qubit state which must be 
designed such as to minimize its presence on the one hand, but to efficiently gather the 
required information on the other hand. 

A powerful readout scheme is based on the phenomenon of dynamical bifurcation, realized 
in form of the Josephson Bifurcation Amplifier {5, 6] and the Cavity Bifurcation Amplifier 
Q . A big J J is placed in parallel to a Copper pair box and driven my an external microwave 
source. Accordingly, close to the first bifurcation threshold determined by the frequency 
and the amplitude of the drive, two stable oscillations appear in the big JJ with thermal 
fluctuations inducing transitions between them. The sensitivity of this process to the shape 
of the Josephson potential is used to retrieve information about the qubit state. However, 
the possibility to tune parameters over wide ranges make these systems interesting on their 
own as devices to study fundamental aspects of driving, nonlinearity and dissipation. In 
particular, the question about the impact of quantum fluctuations on transitions between 
two stable basins of attraction in phase space goes far beyond the standard situation for 
escape over static energy barriers. 

Theoretically, within the relevant range the driven big JJ can be described as a Duffing 
oscillator [7|. This oscillator is particularly important because it represents the simplest 
model to analyze phenomena like bifurcation, period doubling, and dynamical tunneling. 
The classical dynamics of this system is well known and the transition between the two 
stable states, induced by thermal fluctuations, has been investigated in detail js, 9]. With 
lowering temperature quantum mechanical effects appear. It is known that their contri- 
butions are twofold. On the one hand, the transmission probability through the barrier 
becomes finite, on the other hand quantum fluctuations of the environment appear due to 
finite zero point fluctuations in the well adjacent to the barrier and to quantum fluctua- 
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tions in the diffusion coefficients. For escape over static barriers it turns out that for weak 
dissipation the former effect may even exceed the latter one to produce a reduction of the 
escape rate compared to the classical situation [lO] due to a finite reflection of states above 
the barrier energy. Thus, the question arises whether the same is true for transitions over 
dynamical barriers. Moreover, it has to be explored how a consistent semiclassical descrip- 
tion for driven dissipative anharmonic systems must be formulated. To solve both issues is 
the purpose of this paper. 

A powerful procedure for analytical investigations is to describe driven oscillators in a 
frame rotating with a frequency equal to the response frequency of the system (9), ll| . Using 
this approach Dykman and co-workers analyzed the diffusive escape for the Duffing oscilla- 



tor 



12] and for other periodically driven systems 13J. The focus there has been on reservoir 



14j macroscopic quantum tunneling has been 
description of the problem, taking into 



151 ] . In the present work we use a semi- 



induced quantum effects. Complementary, in 
addressed in the deep quantum regime. A numerica 
account multiphoton resonances, has been given in 
classical approach in order to calculate systematically quantum corrections to the classical 
escape rate. We consider both mechanisms, tunneling and bath fluctuations, starting from 
a properly derived master equation. In contrast to 12Ml4|| the impact of finite barrier re- 
flection/tunneling and the whole structure of the dissipative dynamics in the rotating frame 
[l6| are included in this master equation. The corresponding analytical expressions for the 
escape rate apply to the range of weak damping and moderate temperatures. 

The paper is organized as follows: In Sec. [II] we introduce the model and the basic 
notation including the mapping to the rotating frame. This description is extended in 
Sec. IHII to explicitly include also the bath degrees of freedom. This formulation provides the 
basis to derive in Sec. HVI a semiclassical expansion of the master equation the steady state 
distributions of which are used to derive the escape rates in Sec. [Vj In Sec. [TV] this quantum 
diffusion equation is discussed and the quantum corrections to the classical escape rate are 
obtained. 
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II. SYSTEM AND MAPPING ON A ROTATING FRAME 



We consider a system with a weakly anharmonic potential driven by an external time- 
periodic force (Duffing oscillator), namely 

Hs(t) = ^P 2 + \Mulq 2 - irg 4 + Fq cos{u d t) . (f ) 

Accordingly, for the anharmonic coefficient, we assume T(q 2 ) <C Moj\ so that driving is 
almost resonant for 

5u = uq - u d < bj d . (2) 

Classically, when damping is taken into account, two stable oscillations with different ampli- 
tudes and phases appear beyond a bifurcation threshold. The latter one depends on external 
parameters such as driving amplitude F and frequency mismatch 5u. In phase space, these 
two stable states correspond to stable basins of attraction which are separated by an unsta- 
ble domain. Thermal fluctuations may induce transitions between these basins that in turn 
carries information about the global shape of the phase space barrier and the environment. 

Theoretically, the difficulty for a rate description in this kind of system is that the Hamil- 
tonian of the isolated system Hs(t) is time-dependent and, therefore, energy is not conserved. 
However, the dissipative system approaches a steady-state situation such that the reduced 
density matrix takes the form p(t) ~ p(t) cos(cu d t) with an only weakly time- dependent den- 
sity p. For the further analysis it is thus convenient to switch to a rotating frame, given by 
the unitary operator, 

U s (t) = e - iaW , (3) 

where a = ^ |^ (ooaq + j^pj and a) = \J~j^ {^>dq — j^pj are the annihilation and creation 
operators for harmonic oscillators in the system, respectively. The transformation Us{t) 
applied to the coordinate q and momentum p give 

UlqUs = Qcos(uu d t) + — - -Psm(uu d t), (4) 

UJdM 

UlpUs = -u d MQsm(cu d t) + Pcos(u d t) (5) 

with Q and P as new (slowly varying) coordinates. From these equations it is clear that the 
unitary transformation is the equivalent to a rotation of the classical phase space. In the 
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FIG. 1: The Hamiltonian function (from Eq.([6|)) in the rotating frame, for a = 1/27. The energy 
is scaled with uj^Mbujl? . The minimum (m) and the maximum (a) in the figure are the stable 
sates with low and high amplitude, respectively, separated by a marginal state (b). 



rotating frame determined by Us(t) the Hamiltonian reads 



H s = Ul H- i%- U s 



(6) 



= Muj d 5ujL 2 4(g + T _5__!) 




with a length scale L = J d if and the bifurcation parameter 




3F 2 r 



(7) 



a = 



32(u 8uM) 3 ' 



Moreover following a rotating wave approximation fast oscillating terms exp(±mcu c ;t) , n > 1 
are neglected such that a time independent Hamiltonian Hg is obtained. For < a < 4/27 
the rotating frame system exhibits three extrema, whereas the two stable ones correspond 
in the laboratory frame to oscillations with low and high amplitude, respectively. They 
are separated by a phase-space barrier associated with an unstable extremum (see figfT] and 
figj2]). For a < only the low amplitude states exist and for a > 4/27 only the high 
amplitude state remains. 
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FIG. 2: The energies E m (solid), Ef, (dotted) and E a (dashed) of the minimum, the marginal state 
and of the maximum as functions of a. The barrier hight AV = Eb — E m is maximal for a = 
and vanishes for a = 4/27. The energies are scaled with L^uj^Mbto. 

III. SYSTEM AND BATH: MICROSCOPIC DESCRIPTION 

In order to describe quantum dissipation, we explicitly introduce a bath coupled to the 
system, so that the total Hamiltonian is given by 



with a system part as in (pQ), a reservoir part H R and an interaction Hj, i.e. [171, \l 



We switch to the rotating frame with the unitary operator of the composite system 



H = H s + H R + EL i 



(8) 







U(t) EE U S (t)U R (t) 




(10) 



whereas b n and are annihilation and creation operators for harmonic oscillators in the 
bath. In the rotating frame the total Hamiltonian reads 



H = H s + H R + Hj 



(11) 



with Hs as in flH]) and 



uj n m n u d M 



n=l 

P 2 \JL cl 



+ l Q2 + OTJSi^- (12) 

where the new bath parameters read 

Win ~ C n . . 

m n = , U n = LO n - U d , C n = — . (13) 

1 - LOd/Un 2 

With the unitary transformation (flOj) the total Hamiltonian H in the laboratory frame 
OH]) is mapped onto a new Hamiltonian H in the moving frame, composed of the system part 
()6]), a reservoir part (fT2"|) with new parameters and an interaction part (fT2]) . The mapped 
composite system can now be described by techniques applied for undriven escape problems 
with the notable difference though, that the interaction between system and environment 
becomes more complex containing in addition the the conventional position-position cou- 



pling also momentum-momentum contributions. However, as shown in [16|. following the 
standard procedure all bath properties can be captured by a spectral density and by tem- 
perature. Consequently, in the rotating frame we can adapt the approach developed in [lo| 
for transition rates over energy barriers in the energy diffusive limit (weak dissipation). 

One starts from the time evolution of the density matrix of the full compound W(t) 
obeying the Liouville-von Neumann equation ihdW (t) / dt = [H,W(t)} with an initial state 
VK(0). The relevant operator is the reduced density pit) = tr B {W(t)} after eliminating the 
bath degrees of freedom for which a simple equation of motion does in general not exist. 
In case of weak friction and sufficiently fast bath modes, however, progress is made within 
a Born-Markov approximation. One then obtains a master equation which for the present 
case can be cast in the form 16| 



4. fcpl+ ( £w+ ^ + ^ + _^_) w . (14) 

Here operators C xy are defined according to 

/•CO 

£xy[p] = J o dsK' xy (s)[x,[y(-s),p(t)}} 

POO 

+ dsiK^ y (s)[x,{y(-s),p(t)}}, (15) 

with operators y(s) in the Heisenberg representation and {, } denoting the anti-commutator. 
In the rotating frame the force-force correlator functions are defined by 

K — K' 4- iK" 



^(F x (t)F y (0)) p x,y = Q,P, (16) 



where the bath forces read according to ffT2l) 



Vc„x n , F P = V c ra Vn „ . (17) 



Our goal is now to derive from the above master equation f[Mj) a semiclassical equation 
in the energy diffusive limit to determine the leading quantum corrections (order h) to the 
transition rates between the two phase space basins. In case of no external driving the 
analysis in 10] revealed that a conventional position-position interaction between bath and 
system produces in the corresponding diffusion equation only quantum corrections of order 
h 2 . The leading impact of quantum mechanics is thus due to finite transmission through the 
barrier, i.e. tunneling and reflection. Here, however, we will see that while the contribution 
Lpp has a behavior similar to Cqq, the unconventional bath contributions Cqp and Cpq 
yield a supplementary correction of order Ti. 



IV. SEMICLASSICAL MASTER EQUATION 

An energy diffusion operator can be derived starting either from a discrete or a continuous 
spectrum in the well of one the stable domains. The difference is though, that the latter 
procedure conveniently accounts for barrier tunneling near the barrier top, where the density 
of states is basically continuous. It is then convenient to introduce the occupation probability 
of a well state with energy E via 



P{E,t) = "£5{E-E n )p n {t). 



(18) 



n=0 



Here p n is the occupation probability of a well eigenstate with quasi-energy E n and identical 
to the diagonal part of the reduced density matrix in the energy representation. The ex- 
plicit construction of these states follows a type of WKB-recipe as shown below. The time 
evolution equation (THj) can now be represented (see jlO|) as 

R(E')P(E',t) 



P(E,t) 



IF 



-W E 'e 



E,E' 



n(E' 

R(E)P(E,t)~ 
n{E) 



T(E)^P(E,t) 



(19) 
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with u(E) being the frequency of a classical oscillation at energy E and n(E) the density 
of states. Equation ffT9]) captures the incoming probability flux to and outgoing probability 



flux from the state E according to intrawell transition rates 19] 



I r°° 

W E , E , = - 2 \ dm R {{E\H I {i)\E'){E'\H I \E)pt'\ (20) 

II J — oo 

and the reflection probability R(E) from the barrier and the transmission probability 
T(E) = 1 — R(E) through the barrier. 

In the transition rates the system-reservoir coupling appears in the interaction picture 
H T (t) = e i(H s +H R )t/njj ie -i{H s +H R )t/n with p eg _ ^u r j Zr being the equilibrium bath density 

matrix. Note that the unitary transformation ( flOil does not affect the equilibrium density of 
the bath since [Ur, Hr] = 0. The transition rates (1201) can be evaluated explicitly in case of 
the bilinear system-bath coupling as in ffl2l . and one arrives at a golden rule-type of formula 



w Dqq(E - E>) 

W E ,E' = ^2 



\Q qm (E',E)\ 2 + \P qm (E',E)f 



+ DQp{ f 2 E>) 2iS[Q qm (^ E)*P qm (E', E)] (21) 

with Q qm (E', E) = (E'\Q\E), P qm {E',E) = (E'\P\E) /(u d M), and %Z denoting the imagi- 
nary part of Z . The bath correlation functions D xy correspond to bath forces which couple 
to x and y according to f|T6|) and ffTTl) . respectively, 

/oo 
dtK xy (t)e ltE/h . (22) 
-oo 

In accordance with an effectively Markovian dynamics, we consider a purely ohmic envi- 
ronment with a spectral density I{oj) = M^u in the laboratory frame so that the bath 
correlation functions in the moving frame are given by 16| 

Dq Q (E) = 7M[np{E F + E)(E F + E) 

+{n p {E F -E) + l){E F -E)\ : (23) 
D Q p(E) = i^M[-np(E F + E){E F + E)) 

+{np{E F -E) + l){E F -E)\, (24) 

whereas E F = fluid, 7 = 7/4 is the effective friction constant in the rotating frame and 
np(E) = \j{e^ E — 1) is the Bose-Einstein distribution. Interestingly, ( 1231) displays that 
physically two channels in the bath are now open and accessible for emission or absorption 



of quanta, namely, one with energy E F + E and one with energy E F — E 



16] . 



Following the procedure in 10] one arrives with P(E,t) = P(E' ,t)S(E — E')dE' at 
an ^-expansion of ffl9|) in the form 

d 



with Ws(E) = We,e> for E' — E = 5, where 5 is considered as being of order %. The leading 
order terms in the sum above with k = 1 and k = 2 are kept to get the energy diffusion 
equation for finite transmission in the semiclassical limit, i.e., 

d ( _ d 



P(E,t) 



dE \ ^ + 8E^ «W 



Here, the moments of the energy fluctuations read 



P(E,t). (26) 



^=7kiyy Ws(E) - (27) 

To derive semiclassical transition rates from this diffusion equation one must evaluate the 
energy momenta up to corrections of order h and also include consistently transmission and 
reflection coefficients. Bath induced quantum corrections appear due to bath correlations 
D xy that enter the transition rate Ws(E) [see Eq. (1211) ]. Finite tunneling and reflection coef- 
ficients appear explicitly in the diffusion equation (126]) . but must also be taken into account 
in the matrix elements Q qm and P qm that determine the system part in the transition rates 
(pTj) . We note that quantum corrections due to tunneling and reflection are substantial in 
an energy range ~ hu^ around the barrier top, where R(Eb) ~ T(Eb) ~ 1/2 ~ 0(h°). They 
then generate leading corrections in the escape rate that are of order Ti as well. Quantum 
corrections that include combinations of finite transmission and bath induced fluctuations 
are at least of order fi 2 and can be discarded. The strategy we follow in the sequel is thus 
this: in a first step we neglect bath induced corrections and concentrate on the impact of 
tunneling, while in a second step tunneling is neglected and bath fluctuations are accounted 
for. Eventually, the corresponding individual energy diffusion equation derived from (1251) 
are combined to capture both phenomena simultaneously. 
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V. SEMICLASSICAL TRANSITION RATES 



In the energy range close to the barrier top, where tunneling dominates in the temperature 
range considered, the WKB approximation is not applicable because the classical turning 
points to the left and to the right of the barrier are not sufficiently separated. In the 
domain around the marginal point, however, the Hamiltonian Hs can be approximated by 
an inverted harmonic oscillator with barrier frequency [see Appendix (1B3I) ]. and one 
may exploit that the corresponding Schrodinger equation can be solved exactly. The proper 
eigenfunctions are then matched asymptotically (sufficiently away from the barrier top) onto 
WKB wave functions to determine phases and amplitudes of the latter ones. This way, one 
obtains 

(E\Q) = \\(E\Q)- + r(E)(E\Q}+] (28) 



2 



with the matrix elements 



(E\Q) ± = N{E). 



(29) 



'\ird P H s {Q,E) 

containing the action So(E,Q) = fS P(Q',E)dQ' of an orbit starting at the turning point 
Qx and running in time t towards Q with momentum P(Q, E) at energy E. The complex- 
valued reflection amplitude r(E) of a parabolic barrier is related to the reflection probability 
R(E) = \r(E)\ 2 . The normalization follows from (E\E') = 6(E - E') as 



(30) 



In case of vanishing transmission, E « Ej, and R — > 1, one recovers the standard WKB 



wave function N(E) — > yl/(hu) so that it is possible to use (129]) for all energies, provided 
the length scale where a parabolic approximation for the barrier applies is much larger than 



the quantum mechanical length scale ^h/Mub- 

With ( 129]) we calculate in the semiclassical limit the transition matrix elements Q qm (E', E) 

for finite transmission through the barrier. According to the restricted interference approx- 

I I 

imation [20| we only keep the diagonal contributions of forward/backward waves to obtain 
Q qm (E', E) = N (E')N (E) Uq^Y +r(E)r(E'rQ^\ , (31) 
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where Q^(E) = \ dQ(E'\Q)-Q{E\Q} + and Q u Q 2 denote the left and the right turning 
points of the periodic orbit with energy E, respectively. The ^-expansion for the matrix 



elements is thus found to read 
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2l| 



Q 



scl 



+ K (5) 
2 v ' Q 



(32) 



where here and in the sequel the prime ' at energy dependent functions denotes the derivative 
with respect to energy and 

hoj(E) [Qm 



K, 



(6) 



hco(E) 



2mS 



dQe -iSt(Q,E)/?i 

Q(t,E)'e- i5t/h 



Q2{E) 

Qi(E) 



Likewise, we calculate the P matrix element 

P qm (E',E) = N(E')N(E) [- (P^y + r(E)r(EyP s 



(S) 
scl 



(33) 
(34) 

(35) 



with P®{E) = 5^ dQ(E>\Q)-p{E\Q)+ and the expansion 



P. 



(5) 



scl 



P 



(S) 



; p 



>(*) 



cl 



where 



d(6) 
cl 



hcu(E) 



P(Q2,E) 



dPe -i6t(p,E)/n _ 



27TCUdMi5 Jp{Q u E) 

The matrix element in the transition probability (l2Tj) can then be rewritten as 



(36) 



(37) 



2 + |Pqm| 2 « N 4 A + SB, 



S [Q qra {E\ E)*P qm (E r , E)\ « N*C + 5D, 



(38) 



with coefficients A,B,C, and D specified in Appendix lAl 

For escape processes near the bifurcation threshold, the energy level spacings of the 
eigenstates of ([6]) are small compared to hu^. Hence, the following approximation of the 
bath correlations f )23|) and (124"]) applies 



D QQ {E) = ^M{n-E) + 0{h 2 



(39) 



where 



k = h<jJdCoth(hf3uj ( i/2) , 
12 



(40) 



and 

D QP (E) « z ^Mhu d ^ + 0(fi 2 ). (41) 

It is important to note that the lowest order in the ^-expansion of Dqq is of order h°, while 
that of Dqp is of order %. Consequently, as we shall see later, the -Dgp-term in ( 1211 gives 
no contribution to the classical energy diffusion equation [see Eq. (l4"6]) ]. 

Now, using (!38|) . (I39I) . and (HTj) the ^-expansion of the transition probability (j2T|) takes 
the form 

W 5 = ^M^/N^Ak + |^-iv 4 i + kB - 2 ~^^j ■ ( 42 ) 

Close to the energy minimum of the stable domains, the energy level spacing 5 is approxi- 
mately hu m , with i>j m being the minimum local frequency [see App. |B]. Close to the barrier 
top the energy level spacings vanish so that indeed 6 ~ h and f )42|) is a systematic semiclas- 
sical expansion. 



VI. SEMICLASSICAL ESCAPE RATES 

Following the discussion at the end of Sec. IIVI we start in this Section to consider the 
influence of a finite barrier reflection/transmission in presence of a classical reservoir and 
then proceed to analyze the impact of bath induced quantum fluctuations for classical re- 
flection/transmission. Both mechanisms are eventually combined in the preceeding Section. 



A. Finite transmission 



With the transition probabilities at hand, the semiclassical diffusion equation follows 



from ( 12tj|) and ( 12T|) with the semiclassical density of states n(E) 



P(E,t) 



d 



k d 



— C(E) A(E) 1 1 + -— ] R(E) - T(E) 



where 



and 



C(E) = 2 



2dE r 

1 + R(E) 2 
[l + R(E)] 2 

dQ 



hu){E) 

u(E) 
2tt 



as 



P(E,t) 



^) = M,fdQ d -§ + J ^JdP^. 



(43) 

(44) 
(45) 
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The second term in the generalized action (145!) stems from the P matrix element |(n|P|m)| 2 
in (|2T|) . We emphasize that no terms originating from the mixed matrix elements 
(n\P\m) (n\Q\m) appear in (12 ip . 

For vanishing transmission (R — 1, T — 0) one recovers from (14"51) the classical diffusion 
operator 

which (T461) looks like a classical Kramers equation [22] with an effective temperature k. 
Further, A(E) corresponds to an energy relaxation coefficient that takes into account the 
position-position and the momentum-momentum interaction (fT2l) between system and bath. 
As shown in [3] the bath correlation functions Dqq and Dqp are associated with two 
different effective temperatures due to the fact that a detailed balance condition is not 
obeyed in the moving frame. However, since in (1461) [and in (T4"5|) ] the bath correlation 
function Dqp does not play any role, in this regime, it is possible to define the unique 
effective temperature as /cb^cE = k/2. 

Now, the escape rate is determined by the stationary nonequilibrium distribution P s t{E) 
to (T46l) . which is associated with a finite flux across the barrier and obeys the boundary 



conditions P st = for E > E^ and P s t{E) to approach a Boltzmann distribution in the well 
region. Accordingly, one obtains for high barriers 2Vb/n ^> 1 the classical Kramers result 

r cl = ^A(E b ) e _ 2Vb/K 

KIT 

with the well frequency u m (IB2j) . 

The escape rate in the quantum regime including tunneling but no quantum fluctuations 
in the reservoir can now be evaluated also from (T43l) . This diffusion equation formally looks 
like the one already considered in [h| for undriven systems so that we can use the same 
methods to solve it. In the energy range close the barrier top, where tunneling dominates for 
the temperature considered here, the approximate Hamiltonian (1B3|) leads to the parabolic 
transmission and reflection probabilities 

l + exp( L-Jdj 

R = 7"^ f2n(E-V b )\ ■ ( 49 ) 
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(50) 



The quantum partition function in the harmonic well region is given by 

K 00 I/ 2 

Z = ft — 

with Matsubara frequencies v n = nun/h. For vanishing friction, ( 150]) reduces to the known 
result Zqo — l/[2sinh(a; m 7i//c)]. The escape rate follows again from a quasi-stationary 
nonequilibrium state, this time from the quasi-stationary energy distribution P s t(E) of (1431) . 
given by 



scl 



dEn{E)T{E)P st {E). 



This way one gains 



with the coefficients 



sinh(w m ft/fi;) 

1 scl — 7 7~j \ -D J- cb 



5 



1 2Fl 


"l 
2 2 


a, 2 


| + a, 


1 


-0, 


4" 

9. 




2 ~2 


- a, | + 


f + a, 


1 


+ 0, 


4" 
9. 



\ 



27A(E b )(l -^) 2 /k + 36^ 2 



(51) 
(52) 

(53) 
(54) 



8^A(E b )/K 

and the abbreviation 9 = uj^h/ (ttk). The first factor in this rate expression captures quantum 
effects (zero-point fluctuations) in the well distribution, while the second one, \B\ describes 
the impact of finite barrier transmission close to the top. The latter one can actually prevail 
and lead to a reduction of the escape rate compared to the classical situation due to a finite 
reflection from the barrier also for energies E > Eb (fig. [3]). For a fixed 7, the expansion of 
(152"]) for high temperatures is [lol ] 



scl 



r d (1 - he) 



(55) 



with b\ = 1.04 originating merely from the expansion of B. The well partion function leads 
to corrections of higher order in h. 



B. Bath induced fluctuations 



In this section we calculate the impact of the friction terms Cqp and Cpq in ( fl~4j) on the 
energy diffusive decay. According to the above strategy, we assume here to have classical 
transmission and reflection probability, and calculate from (|25|) the first order h correction 
to the classical diffusion equation (|46|) . 
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FIG. 3: Escape rate (|55[) normalized to the classical rate as a function of the bifurcation parameter a 
for hu} d /3/(2ir) = 0.01 (solid), hu d fi/{2it) = 0.05 (dashed) and hu d p/(2-ir) = 0.1 (dotted). For all the 
lines, we use j3L 2 uj d M8u: = 40, 5u = O.Iwq and the dimensionless friction constant f3jML 2 doj = 0.1. 



For R = 1,T = the matrix elements \Qqm\ and |-P qm | are symmetric with respect to 
5, and therefore the terms in the first line of f }2Tj) do not give contributions of order h to 



the diffusion equation 23|, |24j. The only relevant contributions of order % result from the 



Ss[Q qm (E', E)*P qm (E', E)} term. To calculate it, we must take into account also the next 
order term in the expansion of the respective bath correlation function, namely 

E 



Dqp(E) « ijM 



fuo d - — r + ah\E\ 
\E\ 



(56) 



with a 



cosh(hui,il3) — 1 



. We recall that energy level spacings are considered as proportional 



to h. Accordingly, from ( |25|) we obtain the energy diffusion equation, 



P(E,t) 



CO 



-P(E,t). (57) 

Z7T 



where 



A (1) (£) 



IMit^ nA > ^ 



d5 5 2 Z [Q* cl P cl + Q^P*] . 



(58) 



h jo 

For an explicit evaluation of (158]) it is convenient to return to a discrete representation by 
replacing the energy difference 6 with Mu and / dS with J^i^ so that 



(59) 
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2/27 

a 



Am 



FIG. 4: Escape rate (|60p normalized to the classical rate as a function of a for TwJdP / (2~k) = 0.01 
(solid), htddP/CZn) = 0.05 (dashed) and hu)df3/(27r) = 0.1 (dotted). For all the lines, we use 
Su = 0.lu and f3L 2 uj d M5uo = 40. 

This expression is correct for low energies, where the spectrum in the wells is discrete, and 
approximates ( 1581) very accurately for energies near the barrier top. 

In order to reveal the effects of the Cqp and Cpq terms, we calculate the rate of escape 



from (loTlh Following the standard procedure 22| one finds 

r sc i = r cl e^ 2 , (60) 

where 

u - ^_( aK + Ud) / d E , (61) 



CO d K JE m A(E) 

with 6p = hu>d/{nTr). E m and Eb are the energies of the points (a) and (b), respectively, in 
fig. [TJ The integral in (|6"T|) is proportional to the barrier height meaning that ( 16"T|) is of the 
order of Vb(3. It thus gives a significant contribution to the escape rate as depicted in fig. 0] 
for various temperatures. 



VII. DISCUSSION 



In the two previous sections, we have analyzed separately the impact of the two dominant 
quantum effects on the escape rate including contributions of order h. Since corrections due 
to the combination of the two effects in the transition probabilities are at least of order h 2 , 
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a full semiclassical diffusion equation up to order % is obtained by simply adding the two 
results ([5711 and f )43|) . Hence, we get 



P(E,t) 



A ( AC fl - ft 2o A« + Acf Afl + fiWiA « A ) _ r ] « ^ t) . (62) 



The leading order quantum corrections to the escape rate are then found as 



scl 



cl 



l + 0j 



-6a 



(63) 



The first correction is negligible when uj^ <C u,i, i.e. when a approaches the boundaries of 
the bifurcation range (a — > and a — >■ 4/27). Interestingly, the two types of quantum 
fluctuations have opposite effects on the rate expression: while a finite reflection leads for 
energies above the barrier top to a suppression of the escape probability (fig. (j3J)), bath 
induced fluctuations produce an increase (fig. (jSJ)), which typically prevails. The conclusion 
is thus that in the semiclassical regime finite tunneling through the phase-space barrier does 
not play an important role, in contrast to quantum fluctuation induced by the reservoir in 
the moving frame. We recall, that the opposite is true for energy diffusive escape processes 
over static barriers where tunneling leads to a reduction of the rate [10]. 
It is appropriate to remark that the above result is valid for 



(64) 



in order to guarantee the existence of a steady state distribution of a quasi-continuum of 
thermally smeared states on the one hand and to restrict tunneling to energies close to the 
barrier top on the other hand. Equivalently, the range of validity of the rate expression 
is determined by those values of a which are sufficiently smaller than a = 4/27 (where 
AV ->■ 0) [see fig. ©]. 



We have also assumed a weak dissipation compared to the retardation scale H/k of the 
reservoir and to the time scale of the motion 1/utb, l/u) m . However, for a = the barrier 
height stays finite, while goes to zero like a 1//4 and the effective mass M& (see Eq.( 1B3l) ) 
tends to infinity like a _1//2 . The growth of the mass is equivalent to an increase of the 
damping which is mostly clearly taken into account when one introduces a rescaled damping 
constant 7b = 7a -1 / 2 . Hence, for decreasing a, the effective friction 7b/^b grows as well 
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FIG. 5: The dimensionless barrier height AV/(L 2 u)dM5(jj) (solid line), the well frequency u) m /u)^ 
(dotted line) and the barrier frequency ojf, / ui^ (dashed line) as functions of the bifurcation parameter 
a for Soj/iOd = 0.1. 




FIG. 6: Typical behavior of 7b- For a sufficiently bigger than zero 7b is smaller than lo^ and k/%. 
When a goes to zero, 7b grows and achieves a region where it is bigger than uy, but still smaller 
than k/%. For even smaller a, 7b is bigger than Wb and k/%. 

[see fig. (jSJ)] and, therefore, for sufficiently small a the motion near the barrier becomes 
overdamped. 

In the underdamped regime higher order corrections in the friction appear on the one 
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hand through Ub — > Xb in the factor B, where the Grote-Hynes frequency Ab 25| is given by 

Ab = #+^-|. (65) 



On the other hand, the friction dependence of the partition function must be taken into 
account [10j so that 



scl 



n 



l + 9i 



(66) 



n=l "n-Wb + V nl 

with Matsubara frequencies u n = 2irn/%f3. The motion in the well remains always in the 
regime w m /7 < 1. 

To capture the turnover from weak to strong dissipation, one follows the standard Pollak- 



Grabert-Hanggi approach 26|, \2j\ and introduces normal-mode coordinates in the parabolic 



range around the barrier top, where the total system is separable. Then one studies the 
dynamics of the unstable normal mode in presence of the coupling to the stable ones due to 
the potential anharmonicity. For weak dissipation one recovers fl52l) up to leading order. 
With further decreasing a the motion near the barrier becomes overdamped 7b > oJb- A 



classical description applies as long as Ti^b < K [light gray region in fig. 
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a classical Smoluchowski domain 



29| leading to 



30|. Eventually, for very small a, friction becomes so strong 



that ^7b > k [dark gray region in fig. (|Sj)] meaning that the classical Smoluchowski regime 



turns into the quantum Smoluchowski range 
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291 ] . The conclusion of this analysis is that 



almost everywhere in the bifurcation parameter range < a < 4/27 the rate expressions 
( 1631) for very weak friction and its extension to somewhat larger fiction f l66l) are valid. 



VIII. CONCLUSIONS 



In this paper we analyzed the impact of quantum fluctuations on the escape process in case 
of a dynamical barrier and in presence of a dissipative environment. In the energy diffusive 
domain of weak friction and higher temperatures a semiclassical procedure allowed to derive 
effective diffusion equations including leading order quantum effects. It turns out that there 
are two dominant mechanism for these effects to appear, namely, finite transmission through 
the barrier and reservoir induced quantum fluctuations in the moving frame. The latter 
ones dominate by far the deviations to the classical escape rate so that an enhanced escape 
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probability could experimentally be related to the position-momentum coupling terms that 
appear in the rotating frame description of the reservoir resting in the laboratory frame. 

Interestingly, when the bifurcation parameter tends to zero, the strongly underdamped 
dynamics turns into an overdamped motion around the barrier top with a friction strength 
that may even exceed the thermal energy scale. This quantum Smoluchowski domain that so 
far has only been studied for escape over energy barriers (see e.g. 
in a future publication. 



28 



291 ]) will be addressed 
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Appendix A: Coefficients for transition matrix elements 



Here we collect the coefficients appearing in the ^-expansion of \Q 



qm | 



3[Q* P qm 2] in Q3SD. One has 



.4 



+ 



P. 



(5) 



(R 2 + l)+R( Q [ f + Q 



)(sy 



-P. 



(<5) 2 



P. 



B = -(N 4 A)' + N 4 



Q 



cl 



p 



[sy 



' (Q 



cl 



p 



cl 

(<5) 2 



cl 



(R 2 + 1) (Q W + Q^k$) + 2R (g«* « + Q(f^ 



#Q$>P$>* + R{Q®P<1> 



C = N 4 S [(Q 

D = \ {N A C) ' + N 4 [rr*'Q% r P^* -rVQ?P ( 



n {sy p (sy 

Vcl ^cl 



(S) 
cl 



| Pqm | 2 an d 



(Al) 



(A2) 



Appendix B 



Close to its minimum (Q = Q m (a), P = 0), the Hamiltonian ([I]) can be approximated by 

p2 



H, 



(m) 



eS 



2M n 



where the effective mass is determined by 

M" 1 



d 2 H s 



dP 2 



&<* ( Q 2 m 
Mu d { L 2 



(Bl) 
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and the effective potential is V^\Q) = ^MmCo^Q 2 with the frequency 



u m = Sou, 



1 - 3Q r 2 



)2 

' rn 



(B2) 



In the same way it is possible to approximate the system Hamiltonian close to the saddle 
point (Q = Q b (a),P = 0), by 

p2 ] _ 



H. 



(6) 



(B3) 



with M 6 -1 = 5u/(Mu d ) (1 - (Q b /£) 2 ) and 



1 - 3(Q 6 /L) ; 
\ 1 - (W^) 2 



(B4) 
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